Correlations and fluctuations of a confined electron gas 
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Abstract 

J^ ■ The grand potential Q. and the response R = —dil/dx of a phase-coherent 

confined noninteracting electron gas depend sensitively on chemical potential 
C^ ■ fj. or external parameter x. We compute their autocorrelation as a function 

. of /i, 2; and temperature. The result is related to the short-time dynamics 

Q ■ of the corresponding classical system, implying in general the absence of a 

o : 

universal regime. Chaotic, diffusive and integrable motions are investigated, 



^ ■ and illustrated numerically. The autocorrelation of the persistent current of 

a disordered mesoscopic ring is also computed. 
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A major advance in mesoscopic physics has been the observation of single-electron in- 
terference or more generally of quantum mechanical effects in mesoscopic objects. Among 
these we can mention the supershell structure in atomic clusters, the persistent currents in 
mesoscopic rings, or the quantization of the longitudinal conductance in microj unctions [|I|. 

In all these experiments the spatial confinement of the electron gas plays a significant role. 
When a parameter controlling the system is varied - like the chemical potential, a magnetic 
field, or the shape of the confining potential - the response of the gas shows irregular 
fluctuations as a function of the parameter. While the average value of the variation is 
related to a macroscopic property of the system (like its volume), the wild fluctuations are 
the imprint of quantum mechanical interference. 

Neglecting the interactions between the electrons, our purpose is to provide a general 
description of the fluctuating properties of the many-particle ground state electron gas when 
different parameters - like temperature T, chemical potential /z, or external parameter x - 
are varied. Taking as reference parameters (xq, /xq, Tq), the system is described in the grand- 
canonical ensemble by the grand potential 

fi(xo, /io, To) = - j dE X{E, xo) f{E, /io, Tq) , (1) 

with J\f{E,xo) = QsY^k^iE — Ek{xo)) the counting function of the single particle states 
Ek{xo), 9 the Heavyside step function, f{E,fio,TQ) = 1/[1 + exp{{E — /io)/A;To)] the Fermi 
distribution and (7^ = 2 for spin degeneracy. 

The reaction of the system to a variation of the external parameter is characterized by 
the response function 

R{xo,/xo,To) = - dn/dxo\^^j,^ . 

The physical interpretation of R depends on the geometry of the sample and the nature of 
the parameter x. It is a magnetic moment if x is a magnetic field, a persistent current when 
the sample geometry is annular and x is an Aharonov-Bohm flux, or a force when x controls 
the shape of the confining potential (like when deforming a metallic cluster for example). 



We characterize the fluctuating properties of the electron gas by computing the autocor- 
relation function 

Co(x,/x,r) = (f](q_)f](q+)) , (2) 

with q-t = (xo ± x/2,fiQ ± yu/2,To ± T/2). An analogous definition holds for Cr{x, fi,T). 
We denote the fiuctuating part by an upper tilde symbol (see below). The brackets denote 
some energy and parameter averages. 

Parametric correlations like Cq or Cr have been studied in the past for chaotic and 
disordered systems for "individual" (single-particle) levels, rather than for a many-body 
ground state. The single-particle analog of Cr{x, 0, 0) is the velocity- velocity correlation 0. 
Although for the latter no global expression is known, its tail x -^ oo has been shown to be 
universal 10-0. In contrast to this, we will see below that for a Fermi gas the parametric 
correlators are explicitly computable and are expressed semiclassically, even at zero tem- 
perature, in terms of the classical short time dynamics of the system. Because the latter 
is system-specific, the correlators are in general system-dependent with no universal regime 
whatsoever. This difference with the universal local properties occurs mainly because the 
typical time scale of a perturbed ground state Fermi gas is h/fiQ, which is much smaller than 
the Heisenberg time. 

We express the counting function in Eq.(|l|) as a sum of smooth plus oscillatory terms 
A/" = AT + A/". The smooth part is given by the usual Weyl or Thomas- Fermi approximation. 
Semiclassically, the oscillatory part is A/" = '2hgsJ2p{^p/'^p) sm{Sp/h). The sum is over the 
periodic orbits p with action Sp (including the Maslov index), period Tp = dSp/dE, and 
stability amplitude Ap |^. This sum is inserted in Eq.(|T]) to compute the oscillatory part of 
the grand potential. To leading order in h and for kTo ^ /iq 



n(xo, Aio, ^o) ~ '^h'^gs J2 ^p cos(5'p//i) , (3) 

p 

with Ap = Ap k{tp,To)/t^. The factor K,{Tp,To) = gi^hfa/r,) ' ^^^^^ ^^ = h/nkTo takes 
into account finite temperature effects and acts as an exponential cut off for long periodic 



orbits. All the classical quantities entering this expression are evaluated at (xo,/io)- The 
main dependence on x comes from the oscillatory terms via the dependence of the actions, 
Sp = Sp{xo,fio). Taking this into account, the oscillating part of the response function 
R= -dQ/dxo\^,o,To is 

^(xo, yUo, To) ^ 2hgs J2 ^p Qp sm{Sp/h) , (4) 

p 

with Qp = dSp/dxo\f,o,To- 

Replacing Eq.(^) in (|^) we obtain a double sum over periodic orbits with oscillations 
depending on the difference of their actions. The averages in Eq.(§) are done over some 
range A/i and Ax which contain several oscillations of Q but are on the other hand small 
in comparison to the scale over which the classical functions vary. As mentioned before, the 
main dependence on (x, ^) comes from the actions, which may be linearized with respect to 
these parameters. The average restricts the sum to orbits having approximately the same Qp 
and Tp. Ordering the orbits by their period, taking into account these restrictions and using 
the semiclassical definition of the parametric form factor K{t, x) (i.e., the Fourier transform 
of the two-point correlation function B) we obtain 



27r2 io r^ 



Cn = ^ — Xir, T) cosifir/h) Kir, x) , (5) 



with xi'^iT) = k{t,Tq — T/2) k.{t,Tq + T/2). The autocorrelation of the response function 
follows from Eq.@ by derivation with respect to the external parameter. 

Due to the factor r'^ in the denominator, the integral in (^) is dominated by the short- 
time dynamics of the system. For chaotic systems a naive replacement of the form factor 
i^(r, x) by the random matrix result produces a divergence. In real systems, this divergence 
is avoided by the cut off introduced by the shortest periodic orbit of period r^m (the form 
factor vanishes for r < T„iin)- Then even at To = T = (% = !), the autocorrelation will be 
described with good accuracy by computing in Eq.(|^) only the contribution of the shortest 
periodic orbits, the longer ones representing a small correction. Because this short-time 
dynamics is system-specific, then in general we do not expect any universality for Cq, and 
Cr. 



Eq.(|^) may therefore be restricted to the diagonal approximation p = p' for the form 
factor, Kd{t) = h? Y,p ^ ^(t — Tp), with the result 

C^ = ^ —Xcos{fXT/h){cos{Qx/h))Kn{T) . (6) 

zvr^ JO T* 

Here the average is evaluated over the Qp's taken among all the periodic orbits having a 
period between r and r + dr. The expression for the response function is 

Cr = ^2 —Xcos{ixr/h){Q'cosiQx/h))Kn{T) . (7) 

ZTT^ JO T* 

To leading order, Eq.(P) is also valid for systems with a fixed number of particles n [§,0. 
The form factor is now evaluated a.t E = fi{n,xo) obtained by inverting n = Ar(/2,a;o). On 
the other hand, Eq.(^ is modified according to Q^ ^ Q'^ + r^ {djl/dxo) . The importance 
of this new term depends on the nature of the external parameter. If x is a flux in an 
Aharonov-Bohm geometry, AT is independent of x and dfl/dxQ = 0. On the other hand, 
if X is the deformation of a two dimensional billiard, we find that the contribution of this 
additional term is comparable to the Q^ term if the area of the billiard varies, while it is 
smaller by a factor 1/n when the area is kept fixed. 

Eqs.(^ and ([^) are the basic results of this paper. We now proceed to analyze them for 
different physical situations. First we consider a ballistic motion in a regular cavity (the 
trajectories of the electrons occur on phase-space tori). As an example of regular dynamics 
we have studied a Fermi gas contained in a two dimensional rectangular box of fixed area 
A and sides a and b. We choose x = a/b as the external parameter. Exact expressions for 
Ap, Tp, Sp and Qp are known in this case. We can therefore explicitly evaluate Cq and Cr 
as a function of /i or x. In Fig. 1 we show the numerical results for xq ~ 1.3 compared to 
the theoretical curves using only the first 10 orbits as a function of the rescaled variable 
C — a//^ x/xq . We observe erratic long-range oscillations accurately described by our 
formulas. 

We now turn our attention to chaotic systems. Considering that the form factor is strictly 
zero for r < Tmin, the simplest approximation one can made for these systems is to assume 
a random matrix theory behavior for Kd{t) starting at r = Tmm 



, , I ^ ) T ^ Train , , 

KD{r) = { (8) 

2t/(3 , Tmin <T «Th . 

where (3 = 1(2) for systems with (without) time-reversal symmetry. Because the short-time 
dynamics usually dominates the correlation functions, this oversimplified approximation that 
doesn't consider individually the contribution of other short trajectories is not expected to 
be very precise in general. However, what is interesting about Eq.(|^) is that it leads to 
universal expressions for Cq and C/j, since as we will now see the system-dependent features 
may be re-absorb by a rescaling of the parameters. These universal results then serve as a 
guide (in particular in experimental situations where the shape of the potential is not well 
known), and are useful to obtain simple qualitative estimates of the different characteristic 
scales involved. 

Let's evaluate using (H) the variance of the fluctuations of the grand potential and of the 
response function. Restricting to zero temperature, we obtain from Eq.(^ 

C°=C^(0,0,0) = -|^-|^, (9) 

where Tmin is computed at (xo,/io)- This simple expression relates the typical fluctuations 
to the period of the shortest orbit. For a bilhard, Tmin = (mf^l^oY^'^Lmin, with Lmin the 
length of the orbit. Weyl's law allows to express /iq in terms of the average number of 
particles in the gas n = J\f(fj.Q,XQ). We find Cq ~ n'^^'^, with d the space dimensionality. 
This may be compared to the variance of the sum of m energy levels located in a small 
energy window around /iq, which is universal and proportional to m'^logm |l^. If a finite 



temperature is considered, we find that Cq is damped by a factor 4(rmj„/rc) exp{—Tmin/Tc) 

when Tc« Tmin- 

For the variance of the response function we must compute the variance of the distribution 



of the Qp's. We assume for chaotic systems a Gaussian distribution [^,0. The variance is 
(Q^) = a T. Explicit expressions for a valid for chaotic billiards were given in [|13|. Its value 



depends on energy, with a dependence ~ E^^"^. These results as well as the approximation 
(I) in EqM) lead to 



C°=C^(0,0,0) = -^^. (10) 

It is interesting to compare C^ to the "local" variance of the response of the individual single- 
particle levels Ek located in a small energy window around {xo,^o), (f^) = {{dEk/dxY) = 
2agl/ I3th jT^JTsU {th = hp is the Heisenberg time, conjugate to the mean level spacing). 
We have 

Cl = {l/2n'){rH/T^in){v') . (11) 

Because r^ >> Tmm, the electron gas amplifies the local fluctuations. For billiards the 
amplification factor is proportional to n^'^^^^^'^, while the overall dependence on n is C^^ ~ 

The autocorrelation functions are now evaluated in the approximation (^. Setting /i = 
T = To = in (H) and using {cos{Qx/h)) = exp(— arx^/2/i^) by the Gaussian assumption 
we get 



9ii^^^ = (1 - e/2)e-^''/' + eno, em/A , (12) 

where T{a, z) = J^ f^^^ e^* dt. We have rescaled the parameter according to 



e = yS^V^ x/h = ^Cl/2Cl X , (13) 

which makes, at this level of approximation, Cq/Cq a universal rapidly-decreasing function. 
The typical parameter scale over which Cq decorrelates is 6x = h/y/aTmin- This is consistent 
with the previous assumptions made, e.g., 6x is small in the semiclassical h ^ limit. 
It is interesting to compare 6x to the analogous local scale for a single-particle level 0, 



5xl = l/pV< v"^ >. Since 6x = gsJTH/2l3n'^TminSxL, the scale of the Fermi gas is larger by 
a factor ~ j^{d-i)/2d_ 

Expressed in terms of the same rescaled variable, the autocorrelation function of the 
response is 

^^^%^ = e-'^l' - 3 e^r(0, eV2)/2 . (14) 



Analogous results are obtained when computing the autocorrelation in energy, (7^(0, e,0) 
and Ci?(0, e, 0), were the rescaled parameter is e = fiTmm/ft-- 

In order to test numerically the goodness of these approximations we have considered 



an electron gas in a chaotic cavity. The cavity is the two-dimensional Limagon billiard [jl5 
whose shape is defined by the conformal transformation w = z + h z^ + c e^^z^ . For parameter 
values around 6 = c = 0.2 it was numerically shown that the dynamics is dominated by 
chaotic trajectories [ll5i[T6| . We take x = 5 as the external parameter. Fig. 1 shows the 
numerical results compared to the theoretical curves. The approximation (|^) provides a 
good qualitative description of the main features. The deviations observed may be due to 
contributions of periodic orbits other than the shortest one. 

Analogous results may be obtained for diffusive systems by using the appropriate form 
factor in Eqs.@ and (|^. We concentrate on a particular example, the calculation of the 
autocorrelation of the persistent current of an electron gas in a disordered ring threaded by 
an Aharonov-Bohm fiux 0. 

The current at thermodynamic equilibrium is given by J = — c^. Only the fiuctuating 
part of fi depends on 0. We replace S^^jli -^ Sp/Ti + 27rWp0/0o in Eq.(^ to include the fiux, 
with 00 = hc/e the fiux quantum and Wp the winding number of the periodic orbit p around 
the ring. The current in the ring is 

/(A) = — - — - ^ Ap Wp sin {Sp/Ti + 2T:Wp\) , (15) 

00 p 

where we have introduced A = 0/0o- The autocorrelation of the current C/(A, \r) = (/(Ar + 
A) I{Xr)) is computed in the diagonal approximation taking into account the contribution 
of the primitive orbits and of their time reversal partners. As before, it can be expressed in 
terms of the form factor, which now takes into account the diffusive motion of the particle 
|T^,|T8|. For simplicity we restrict the calculation to a quasi one-dimensional ring and set the 



temperature to zero. At = we get the well known variance of the persistent current C^ 
|T^,[TB|. Normalizing Cj with this variance and averaging both with respect to the reference 



flux Aij, we arrive at 



CijX) _ 1 ^ cos(27rwA) 
C? C(3)i, t.3 ' ^""^ 

where ({s) is the Rieinann zeta function. 

Eq. (|T6|) is a universal function, independent of the specific aspects of the dynamics (i.e., 
independent of the elastic time between collisions or the diffusion constant), but depends on 
the dimensionality of the system. No rescaling of the parameter was needed. This interesting 
fact is due to the geometry of the sample and the physical quantity considered. Orbits with 
zero winding number do not contribute to the current because they have no flux dependence, 
while those with w > 1 have a period t >> r^, the elastic time. As a consequence, all the 
non-universal shortest orbits gave no contribution, and the autocorrelation is a universal 
function (at fixed d). 

In conclusion, we have studied the statistical properties of a phase-coherent noninter- 
acting confined electron gas in response to a parameter. The results are applicable to 
any confined fermionic system. The autocorrelations Cq and Cji were computed explicitly. 
Their behavior is described by the shortest classical periodic orbits, making them system- 
dependent with in general no universal regime. When the form factor is simplified to take 
into account only the shortest orbit, it is possible to remove the dependence on specific 
properties of the system by a rescaling determined by Tmm- The typical fluctuations of 
the gas are also expressed by simple formulas in terms of Tmin- This is different from the 
rescaling of parametric correlations of single-particle levels, determined by th. This latter 
time is of secondary importance in our context, and long orbits represent a small correction 
to our results. In the case of persistent currents in a diffusive Aharonov-Bohm geometry, 
the autocorrelation is independent of the parameters of the system, without any rescaling, 
due to the suppression by the geometry of the short-time orbits. The analytical results have 
been verifled by numerical simulations on chaotic and regular cavities. 

We are especially grateful to H. Bruus, C H. Lewenkopf and E. R. Mucciolo for providing 
the program for the numerical simulation of the Limagon billiard. 
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FIGURES 




FIG. 1. Normalized autocorrelation functions (points) compared to theoretical predictions 
(dashed lines): (a) Cq{£,,0,0)/Cq and (b) Cr{^,0,0)/Cj^ for the (integrable) rectangular billiard; 
(c) Cn(0,e, 0)/(7q and (d) (7^(0,6, 0)/C^ for the chaotic Limacon billiard 
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